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We investigate the electrical switching of charge and spin transport in a topological insulator 
nanoconstriction in a four terminal device. The switch of the edge channels is caused by the coupling 
between edge states which overlap in the constriction and by the tunneling effects at the contacts 
and therefore can be manipulated by tuning the applied voltages on the split-gate or by geometrical 
etching. The switching mechanism can be conveniently studied by electron interferometry involving 
the measurements of the current in different configurations of the side gates, while the applied bias 
from the external leads can be tuned to obtain pure charge or pure spin currents (charge- and spin- 
bias configurations). Relevant signatures of quantum confinement effects, quantum size effects and 
energy gap are evident in the Fabry-Perot physics of the device allowing for a full characterization 
of the charge and spin currents. The proposed electrical switching behavior offers an efficient tool 
to manipulate topological edge state transport in a controllable way. 



I. INTRODUCTION 

The discovery of Topological Insulator systems (TIs), 
both in three and in two dimensions, has recently at- 
tracted enormous attention^. TIs possess an insulat- 
ing bulk gap and metallic edge or surface states, which 
can be distinguished from an ordinary band insulator by 
the existence of topological invariant 3 -. This excep- 
tional property leads to quantum spin Hall (QSH) ef- 
fect which was first proposed for a model graphene sys- 
tem by Kane and Mele 4 -. As a specific example of two 
dimensional QSH system, a HgTe/CdTe quantum well 
(QW) with an inverted band structure has been demon- 
strated, experimentally and theoretically, to have a single 
pair of helical edge states in the QSH bar by the ap- 
pearance of a quantized conductance plateau when the 
Fermi energy lies in the bulk gap 5-7 . Quantized transport 
along the HgTe boundaries can be conveniently explained 
by an edge channel picture: Two states with opposite 
spin orientation propagate along opposite device edges 
in the same direction thus leading to a quantization 8 - of 
the conductance of 2e 2 /h. Due to the spatial separa- 
tion of the spin-states in these systems and to their one- 
dimensional (Id) nature, system geometry and interfer- 
ence phenomena can be conveniently used for spin selec- 
tions and HgTe-based topological insulators appear to be 
promising candidates for spin processing devices. Recent 
proposals of spin-transistors based on two-dimensional 
topological insulators rely on the application of a mag- 
netic field at a pn-junction 9 , Aharonov-Bohm and Fabry- 
Perot interferometers^—, and gating of a single HgTe 
nanoconstrictio n 13 i 14 . 

In this paper we demonstrate how topological edge 
states can be electrically switched in an elongated con- 



striction leading to charge and spin transport with high 
fidelity. The physical mechanisms utilized here involve 
the coupling between topological edge states (TESs) and 
the tunneling effects (including spin-flip tunneling) at the 
extremes of the constriction due to local etching. Both 
the coupling between the TESs and interference phenom- 
ena along the constriction can be controlled by all elec- 
trical gating and the present analysis sheds light on how 
manipulate edge-state transport in TI. 
In the following we introduce a Id effective model. Usu- 
ally the QSH physics in HgTe/CdTe QW was studied 
by an effective 4-band model that depicts the inversion 
crossing of electron and hole band 5 and most of the works 
have been based on the numerical solutions within a 
tight-binding method^ 5 -, while an analytical solution for 
the case of a finite strip geometry recently appeare d 16 ' 18 . 
The transverse finite size effect is relevant because the 
edge states on the two sides can couple together to gen- 
erate a gap in the spectrum even in the clean limit, break- 
ing the edge channels. Since a gap opens in the spectrum, 
differently from the quantum Hall edge states which do 
not couple across the width of the strip, the application 
of electrical gate potentials can be exploited to shift the 
position of the electrochemical potential within the gap, 
thus permitting the switching of spin and charge trans- 
port. The transport properties of QSH systems in pres- 
ence of an extended contact have been also considered 
in the interacting case both to extract information about 
the intensity of the interaction^ 9 , and to investigate the 
limit of extremely narrow constrictions^. 
The organization of the paper is the following. In SecHIl 
we introduce an effective one-dimensional model to de- 
scribe edge states at the surface of a 2D topological insu- 
lator and their coupling along a nanoconstriction. Here 
we also present the model for the four terminal set-up and 



its operational configuration. In Sec lIIII we introduce the 
scattering field approach and the basic formalism used 
to calculate spin and charge current in terms of the scat- 
tering matrix elements. In Sec lIVI we show the results 
concerning charge and spin currents, mainly focussing 
on two specific configurations of the four terminal setup. 
Finally, we discuss our conclusions in Section fVl 

II. MODEL HAMILTONIAN 

We consider a QSH bar with a nanoconstriction of 
transverse dimension W and length L formed by a split 
gate or by a geometrical etching [see FigJI]. Starting from 
the 4-band model of Refs.0-0] one can derive an effec- 
tive Hamiltonian for an infinite strip of constant width 
W. In the region close to the external leads the strip is 
characterized by a wide transverse dimension larger than 
the transverse decay length of the edge state wavefunc- 
tion. Here one can derive an effective Dirac Hamiltonian 
for the edge state -ff e //,t/J. = c ~F h v F&xk x ipx being the 
Pauli matrix) in a single spin sub-block with a velocity vf 
which corresponds to the Fermi velocity and energy offset 
c which is in agreement with the full band structure in 
the vicinity of the band crossing^. For decreasing width 
W, i.e. along the constriction, the edge states at oppo- 
site boundaries start to overlap, leading to a mass like 
gap in the particles spectrum whose size is an exponen- 
tially decreasing function of the width T cx exp(— XW), 
as explained in Appendix [XI while A is determined by the 
secular equation for the eigenvalues and depends on the 
distribution of the wavefunctions in spaced. Addition- 
ally, one could also take into account the leading order 
spin-orbit interaction (SOI) due to bulk or structure in- 
version asymmetry but the overlap of edge states due 
to both of them is negligibly small and is relevant only 
close to the avoided band crossing and for constrictions 
of width of several tens of nanometers (< 100 nm). We 
can thus consider the following Id effective Hamiltonian 
for a Kramers pairs edge states along the quantum well: 

H = H + H sp + H sf + H C1 (1) 

where 1 ': 

H a = -ihv F y^ dx [: ij)\ la .(x)d x il) Ra {x) : 

-:^ Ls {x)d x ^ Ls (x) :], (2) 

and iPr(l)<t represents the right (left) mover electron an- 
nihilation operator with spin a =f, 4-, while : O : stands 
for the normal ordering of the operator O with respect 
to the equilibrium state defined by occupied energy lev- 
els below the Fermi sea. In our description, without loss 
of generality, we assume that spin-| right movers (R, f) 
and spin- J, left movers (L, J,) flow along the top bound- 
ary while the spin- .J, right movers (R,i) and spin-f left 
movers (L, f) flow along the bottom boundary. Along 




FIG. 1: (Color online) Representation of the edge states flow 
at the top and bottom boundaries of a TI quantum well. In 
the middle region, characterized by transverse dimension W 
and length L, a finite overlap among the surface states allows 
the formation of a gap in the quasi-particle spectrum, while, 
close to the leads, states with different helicity are uncoupled. 
The transport properties of the system can be controlled by 
side gates, i.e. Vt and Vb- In particular two gate configura- 
tions are analyzed: (i) Vt = Vf, = V g ; (ii) Vt = —Vt = V g . 
Two bias configurations are considered: (i) Charge-bias de- 
fined as Vi = Vi = V, Vz — Vi — 0; (ii) Spin-bias defined as 
-Vi = V 2 = V, V 3 = Vi = 0. 



the nanoconstriction the edge modes are coupled by con- 
finement effect and the overlap between edge states be- 
longing to different boundaries open a gap in the energy 
spectrum of the Dirac Fermions. At the extremes of the 
constriction (x ~ x\ and x = X2 ) inter-boundary tun- 
neling events may take place and the only terms which 
preserve time-reversal symmetry^ can be distinguished 
in spin-preserving and spin-flipping tunneling described 
by the following Hamiltonians: 




+T sp {x)*ip\ g {x)ip m {x)}, (3) 




+T sf (xr^ ai (x)Mx)}, (4) 

where a = {L,R}, = +1,£l = — 1 is the chirality, 
while r s j( sp )(a;) are the space-dependent tunneling am- 
plitudes: 

^sf(s P )(x) = 2hv F 2J lsf(s P )S(x - x 4 ). (5) 

i£l,2 
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Finally the term H c in ([1) describes the coupling between 
the edges along the constriction: 

Hc = E / dX T ( X )A( X )^La(x) + 

+ipU{x)i>m(x)}, (6) 



<r=t,4. 1 



where r(x) = C(x)T and C(x) is a step- like function tak- 
ing value 1 along the nanoconstriction (i.e. x\ < x < X2) 
and zero elsewhere. Differently from Ref.[l3j], we as- 
sume that the main spin-flipping mechanism in our model 
is caused by the local modification of the spin-orbit 
coupling- 2 - at x = X\p, governed by j s t, while we disre- 
gard the spin-orbit interaction eventually present along 
the constriction Tt. The latter assumption is fully justi- 
fied for not too tight W where Tf/T < A 

In the presence of side gates Vt, V& at top and bottom 
boundaries of the device, an additional term appears in 
the Hamiltonian: 



H g = I dx[eV t {pR-\ + pLi) + 



(7) 



where p aa —: ip^ipaa : denotes the electron density 
with a = L,R and spin a. The presence of such volt- 
ages shifts the edge state momenta in the top and bot- 
tom region between x\ and xi, modifying the electron 
phase in the loop processes induced by the tunneling and 
can give rise to electron interference phenomena rem- 
iniscent of the Fabry-Perot (FP) and Aharonov-Bohm 
(AB) quantum phases, with respective value given by 
irtfipp = e(Vt+Vb)L/HvF andn4>AB — e(Vt — Vb)L/hvF^. 
In the subsequent analysis we consider different bias and 
gate configurations 2 - resulting in four possible opera- 
tional modes of the device. Concerning the bias ap- 
plied to the four terminals (see Fig. [IJ, we define the 
charge-bias configuration (CBC) with Vi = V2 = V, 
V3 = V4 = 0, and the spin-bias configuration (SBC) cor- 
responding to —Vi = V2 = V, V3 = V4 = 0. The defini- 
tion of 'charge' and 'spin' configuration originates from 
the degree of freedom injected through the scattering re- 
gion, depicted in FigfTJ Indeed in configuration CBC the 
amount of spin-up and spin-down electrons injected from 
terminals 1 and 2 is the same, so that only the charge de- 
gree of freedom is injected, and no net spin. In contrast, 
in configuration SBC the lead 1 is negatively biased, de- 
termining a depletion of spin-down electrons with respect 
to the equilibrium situation, supplying a spin degree of 
freedom to the arrival leads (i.e. the leads 3 and 4). The 
same configurations have been discussed in Ref.|ll| in a 
model without coupling between the edges. In interme- 
diate situations, i.e. when |Vi| ^ | V2 j, both charge and 
spin degrees of freedom are involved, but for clarity we 
focus only on CBC and SBS where only pure charge or 
pure spin current can be generate. Moreover, two side 
gates configurations are analyzed: (i) V t = V h = V g 
(4> FP ^ and cj) AB = 0); (ii) V t = -V b = V g {<f> FP = 



and 4>ab 7^ 0), which permit us to analyze electron in- 
terferometric phenomena. 



III. SCATTERING FIELDS APPROACH 

We now formulate a scattering field theory a la 
Biittiker^ able to describe coherent spin and charge 
transport in the system shown in FigfT] 
The charge or spin current operators J c / S in first quan- 
tization are written as follows: 



Jc 



v F er z 

h 

v F -r z 



(8) 



where a z (f z ) stands for the Pauli matrix, I2X2 for the 
identity matrix acting on the Hilbert space given by the 
tensor product \a) ® \a) (a € {R, L}, a E {t,l})- To 
build a scattering field theory, one first defines the scat- 
tering field corresponding to each terminal i = 1, . . . , 4 in 
terms of the incoming (a acr (E)) and outgoing {b aa {E)) 
electron operators, according to: 



*i(x,t) 

$2(x,t) 

* 3 (M) 
* 4 (z,t) 



Et/h 




a Rl (E-x)+b Lt (E-x)\ (9) 
aaf(E;x) +b Li (E\x) 
a Li (E;x) + b m (E\x) 
a L f(E;x) + b Ri {E\x) 



where a atr (E;x) — a a(T (E)\a) ® \a) exp(in a kEx) with 
rju = —r\L = 1 and the wavevector Ue = E/(hv F ) (and 
similarly for 6) . 

The second-quantized current operators in the terminal 
i is defined by = J c / S ^i and is explicitly given by 
e {c,s}): 



i9n 



(10) 



where g c = \e\/h, g s = l/(47r), £ c / s = Tl and ei, 4 = 
— I = —£2.3- In writing Eq. (jlOp we made use of 
the Fourier transform cn(t) = J dE&i(E) exp[— iEt/H], 
while the following correspondence has been made: 
[61,62,63,64]* = [bft, Hi, buf, buxf, [ai,a 2 ,a 3 ,a 4 ] t = 
[o-Rl, flflf; a-Li, Oitl*' 

The expectation value ( ) can be computed making use 
of the scattering relation bj = J^. Sjiai and from quan- 
tum statical average (a^(E)a i (E')) = S tJ S(E - E')f i (E) 1 
being fi(E) the Fermi-Dirac distribution with electro- 
chemical potential pj = p + eVj . After direct computa- 
tion we get: 



dE i^T +1 + (Wio (E)\ 2 } fj (E) 



(11) 
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In the linear response regime fj(E) can be expanded 
around the equilibrium energy e = E — jl and for small 
bias one can express the charge/spin current in terms 



of a generalized conductance tensor 
whose elements are given by: 



J2j g^Vj, 



CP 



de 



(d s m) . 

\ / eq 

(12) 

Eq. (fT2|) provides a description of the linear response the- 
ory of the system in terms of the scattering matrix ele- 
ments. Since we are interested in the quantum regime, 
we shall limit our analysis to the zero temperature case. 

A. Boundary conditions and scattering matrix 



to half gap in the coupling region for the parameters 
here used, while the distance L is rendered dimension- 
less by the substitution L — > EL/(Hvf) = d (notice that 
L w 0.49 fim xd, for v F 0.48 • 10 6 m/s). Finally, we 
take the coupling T = 2 corresponding to 1.3 meV, being 
it an appropriate value for a nanoconstriction of W = 100 
nm (see Appendix The currents are measured in the 
terminals 3 and 4, which are assumed to be grounded, 
i-e. /a = fi = feq, with f eq the Fermi distribution at 
equilibrium, while the J c / S vs d curves are shown on a 
very large range (i.e. d £ [0.2,6]) to better identify the 
characteristic oscillation scales. 



Gates configuration 1: Vt — Vb — V g and bias 
configuration CBC 



The scattering matrix Sij is a four by four unitary 
matrix whose diagonal entries vanish by helicity and 
time-reversal symmetry while all the other entries can 
be explicitly determined as a function of the tunnel- 
ing amplitudes J sp ,j s f by imposing the proper bound- 
ary conditions (BCs) on the wave function a 11 ' 12 . Since 
the wave functions are continuous in the regions x < x%, 
x\ < x < X2 and x > 12, we only have to impose 
the matching conditions where Dirac delta potentials are 
present, i.e. at x = x 1 / 2 . By using the equation of 
motion of the quantum fields (see Appendix [SJ and ex- 
plicitly taking into account the properties of the Dirac 
delta potential under integration, one obtains the follow- 
ing matching conditions: 



MV(xi + 0+) = M*^{ Xl - 0+) 
(i £ {1, 2}) where the matrix M. is given by: 



(13) 



/ 



M 



1 



V 



ilsf ilsp 

1 ij sp 

1 i~f sf 

-ilsp hsf 1 



(14) 



In the limit of vanishing 7^ and 7 sp , M becomes the 
identity matrix I4X4 and thus the BCs simply require 
the continuity of wavefunctions in x = Xi. The BCs in 
(Ti~3l provide 8 equations from which the scattering ma- 
trix elements can be numerically determined. Once the 
scattering matrix is known, the charge and spin currents 
can be computed by using Eq. (1111) . 



IV. RESULTS 

In the following analysis we study the charge and spin 
currents induced through the system as the effect of the 
applied bias Vi, i £ {1,...,4}. We will work in dimen- 
sionless units. In particular, the charge (spin) currents 
are expressed in unit of Ve 2 /h [eV/(47r)], the energy is 
measured in unit of S = 0.65 meV which corresponds 



In this configuration the momentum of the edge states 
along the nanoconstriction is shifted by the Fabry-Perot 
phase and the energy spectrum, which can be obtained 
by a straightforward diagonalization of the Hamiltonian 
Hq + H c + H g , is characterized by two branches: 



E± = eV g ± yJ{hv F k) 2 +r 2 , 



(15) 



originating from the coupling T of left- and right-movers 
in each spin channel. When the Fermi energy is located 
within the gap, for sufficient long L, the particles trans- 
port is strongly suppressed, while applying the gate V g 
the edge channels transport can be switched. In the fol- 
lowing we consider the particles transport through the 
conduction band, i.e. the one specified by the " + " sign 
in Eq. ([T5]) . 

In FigJ3] we present the results for the CBC where a 



2.0 



1.5 



J c 1.0 



0.5 



0.0 




-6 



-5 



-4 



-3 



-2 



e v„ 



FIG. 2: (Color online) Charge currents J c vs eV g in gate 
configuration 1 computed by fixing the model parameters as 
follows: d = 2, r = 2, j s f — 0.1. The parameter "/ 3p varies 
from top to bottom j sp £ {0.05,0.1,0.2,0.25}. The above 
curves, for a fixed value of ' 
7s/- 



do not depend on the value of 



pure charge current J c is generated. By fixing the sys- 
tem parameters as d — 2, T = 2, j s f = 0.1, we plot- 
ted a set of curves of J c vs eV g at varying 7 sp in the 



5 



set {0.05,0.1,0.2,0.25}. When eV g = the conductance 
is strongly suppressed, while for a finite gate potential 
a non-vanishing current can be generated. The J c vs 
eV g has an oscillatory behavior associated to the Fabry- 
Perot-like resonance in the constriction. The maxima of 
J c are determined by the constructive interference con- 
dition along the nanoconstriction: 

2kd = 27rn - 2<j> s , (16) 

where <f> s is a scattering phase not known a priori and 
which depends on the transparency of the barriers, i.e. 
on the tunneling amplitude "f spi j s f at xi.2, while the 
particles momentum k is determined by Eq. (ITBl . The 
values of the bias maximizing J c are instead given by 
(eVg) n = — ^JY 2 + {(mi — cj) s )/d] 2 in correspondence of 
the Fermi level Ep — 0. The dependence of J c on the 
scattering phase </> s is evident in the curves of Fig |2] where 
7s P is varied from lower to higher values. In particular, 
a linear shift of the interference maxima (minima) ac- 
companied by an higher amplitude modulation (at high 
values of -f sp ) is observed. The shifting of maxima of J c at 
varying 7 sp is not observed in the case of Dirac Fermions 
with linear dispersion relation (E cx k) and thus the ob- 
served s -dependent shift is a peculiar feature of gapped 
Dirac particles. In fact, in this case the non-linear en- 
ergy dispersion leads to scattering properties analogous 
to those of massive Schrodinger particles. An experimen- 
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FIG. 3: (Color online) Charge currents J c vs d in gate configu- 
ration 1 computed by fixing the model parameters as follows: 
eV g = —2.45, r = 2, y s f = 0.1. The parameter 7 sp varies 
from top to bottom in the range 7 sp € {0.05, 0.1, 0.2, 0.25}. 

tal study of the oscillating behavior of the J c vs eV g could 
provides relevant information on the coupling energy T 
and on the scattering phase <p s . 

In Fig. |3]we study the charge current as a function of the 
dimensionless length d of the constriction by fixing the 
remaining parameters as follows: eV g = —2.45, T = 2. 
7 s/ = 0.1, while 7 sp takes the values {0.05, 0.1, 0.2, 0.25}. 
Apart from the shifting of the interference maxima, a 
general periodic behavior is observed. The space period 
t ( i depends on the Fermi energy Ep, on the applied volt- 



age V g and on the coupling T and can easily be deter- 
mined by (H5]) and ([TBI) : 

Let us note that does not depend on 7 sp nor on the 
scattering phase <p s which is instead involved in the inter- 
ference condition (|16[) , providing a shift of the maxima of 
J c . Furthermore, the expression of r<j can be experimen- 
tally used to determine the average value of the coupling 
r when more devices with different channel length d are 
at disposal. 

B. Gates configuration 1: Vt = H = V g and bias 
configuration SBC 

The device under consideration (FigllJ can also 
work in a different bias configuration allowing for the 
generation of pure spin current (SBC). In the following 
we present results in this configuration. In Fig|4]we show 
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FIG. 4: (Color online) Spin currents J s vs eV g in gate configu- 
ration 1 computed by fixing the model parameters as follows: 
d = 2, r = 2, 7s/ = 0.1. The parameter 7 sp varies from 
bottom to top curve as y sp g {0.05, 0.1, 0.2, 0.25}. 

the spin current J s as a function of eV g by setting the 
model parameters as: d = 2, V = 2, j s f — 0.1, while the 
different curves correspond to -f sp € {0.05,0.1,0.2,0.25}. 
The oscillatory behavior as a function of V g and the peri- 
odicity of the curves follows strictly the one described for 
the CBC by Eq. (fl~6"| . However, differently from that case 
a lowering of the maxima of the spin currents is observed 
by decreasing 7 sp . This effect is also evident in FigJSJ 
where J s vs d curves are shown by fixing eV g = —2.45, 
while maintaining the other parameters as given in FigU) 
The lowering of the maxima as a function of 7 sp can 
be understood by observing that the expectation value 

of the spin current is given by VF{h/2)\J \ — (T/E) 2 
which depends on the scattering phase 4> s through 
the energy evaluated at eV g . Differently from the 
charge current analyzed in CBC, the spin current J s 
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FIG. 5: (Color online) Spin currents J s vs d in gate configu- 
ration 1 computed by fixing the model parameters as follows: 



eV g 



-2.45, r = 2, 7 S / = 0.1. The parameter -y ap varies 



from bottom to top curve as -y 3 p £ {0.05, 0.1, 0.2, 0.25}. 



FIG. 6: (Color online) Spin currents J s vs d in gate configu- 
ration 1 computed by fixing the model parameters as follows: 



eV g = 



-2.45, r = 2, lap = 0.1. 
from bottom to top curve as j a f 
also APPENDIX [5] 



The parameter -f a f varies 
<= {0.05,0.1,0.2,0.25}. See 



depends on the value of j s f, the spin- flipping tunneling 
amplitude. This dependence is evident in Figj6l where 
the spin current J s is analyzed as a function of the 
constriction length d by fixing the model parameters 
as eV g — —2.45, r = 2, j sp = 0.1, while the different 
curves are computed for j s f £ {0.05,0.1,0.2,0.25}. As 
already seen in FigJSJ an oscillating behavior of the J a 
vs d curves is observed, the period of the oscillation 
being described by Eq. (|17|) . However, an important 
feature in the curves in FigjB] is that by increasing j s f 
above a threshold value (j s f 0.2) the spin current 
changes sign (see APPENDIX |B|. Since 7 S / can be 
locally controlled by geometrical etching, the change 
of sign of the spin current can be implemented as a 
switch for spintronics purposes or to characterize the 
constriction. The above method is expected to be quite 
robust against decoherence phenomena which can reduce 
the Fabry-Perot oscillations amplitude while unaffecting 
the mean value of the current. 




FIG. 7: (Color online) Spin currents J 3 vs eV a in gate configu- 
ration 1 computed by fixing the model parameters as follows: 
d = 2, r = 2, 7s P = 0.1. The parameter 7 S / varies from 
bottom to top curve as j s f £ {0.01, 0.05, 0.1, 0.15}. 



Gates configuration 2: Vt = — V& 
configuration CBC 



In FigJT] we analyze the Fabry-Perot oscillations of 
the spin current J s as a function of eV g by fixing the 
model parameters as d — 2, T — 2, j sp = 0.1, the pa- 
rameter 7 s j being fixed from bottom to top curve as 
jsf e {0.01,0.05,0.1,0.15}. As evident from the figure, 
a decreasing of 7 S / produces an increasing of the spin 
polarized current generated in the system, the maximum 
value being related to the number M of modes involved 
in the transport (i.e. Af = 2). 

As a final remark, we observe explicitly that the charge 
current in FigJ5] (CBC) and the spin current in FigJT] 
(SBC) can be controlled by using V g . Indeed, starting 
from eV g = —2 (off-state), a sudden activation of the 
transport along the edge can be induced by tuning eV g 
to the value —2.4 (on-state). 



In this configuration the side gates energy eV g removes 
the spin degeneracy leading to the following energy spec- 
trum: 



E, 



/3^J(hv F k + rieV g ) 2 +T 2 , (18) 



where (3 = ±1 indicates the conduction (" + ") and the 
valence (" — ") band, while ?/ = ±1 describes the mo- 
mentum shift of particles. The spectrum (|18j) looks very 
similar to the one originated by the presence of spin- 
orbit coupling along the constriction. In particular, the 
eigenstate corresponding to r\ = 1 acquires the phase fac- 
tor exp[— ieV g /(hvp)) compared to the case with V g = 0, 
while that for rj = — 1 the factor exp[ieV g /(hvp )]. In the 
following discussion, we fix the Fermi energy Ep — 2.45 



and r = 2 to obtain a non-vanishing particles transport 
through the conduction band (f3 = +1). The value of Ep 
is a tunable quantity^ and can be controlled by using 
an additional back-gate acting below the whole nanos- 
tructure, the more interesting regime being the one with 
Ep just above the gap. In Fig|8] we study the charge 
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FIG. 8: (Color online) Charge currents J c vs eV g in gate 
configuration 2 computed by fixing the model parameters as 
follows: d = 2, r = 2, -y sf = 0.1, E F = 2.45. The parameter 
7s P varies as j sp G {0.05, 0.1, 0.2, 0.25}. The above curves, for 
a fixed value of 7„ p , do not depend on the value of y s f. 



current J c as a function of the applied gate eV g by fix- 
ing the model parameters as follows: d = 2, T = 2, 
7s/ = 0.1, E F = 2.45, while -y sp G {0.05,0.1,0.2,0.25}. 
A distinctive feature is that V g does not induce a charge 
current modulation as the one observed in FigJ3] This 
phenomenon can be explained observing that the inter- 
ference effects responsible for a charge current modula- 
tion are related to the Fabry-Perot phase <j)pp oc Vt + V& 
which is zero in this configuration. The analysis of J c as 
a function of the constriction length d is shown in Fig[§] 
setting the model parameters as follows: eV g = —2.4, 
r = 2, 7s/ = 0.1, Ep = 2.45, 7sp G {0.05,0.1,0.2,0.25}. 
The J c vs d curves show an oscillating behavior whose 
space period can be deduced by the values of k obtained 
from Eq. fT5]) fixing the energy at Ep. In this way we 
obtain two different space periods: 



y/E F - T 2 ± eV g 



(19) 



.(±) 



which define two characteristic frequencies u>± — 2tt /t\ 
that determine the harmonic content of the charge and 
spin current curves. In particular, as in the familiar case 
of superposition of waves with different frequencies, one 
expects to see an oscillating function with frequencies 
f2i/2 = (w+ ± w_)/2 and whose corresponding space pe- 
riods are: 



r (2) _ 



eV a 



(20) 



2vr 



In FigHthe period /2 
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FIG. 9: (Color online) Charge currents J c vs d in gate config- 
uration 2 computed by fixing the model parameters as follows: 
eV g = -2.4, T = 2, 7 s/ = 0.1, E F = 2.45. The parameter -y sp 
varies as 7 sp G {0.05,0.1,0.2,0.25}. 



while the period ~ 1.31 is not observable being near 

(2) 



equal to a subharmonics of /2 



D. Gates configuration 2: Vt = — Vb 
configuration SBC 

When the system is driven in the SBC, pure spin cur- 
rents are generated. In the following we focus on this 
configuration. 

In Fig 1 101 we present the spin current J s as a function of 
eV g by setting the model parameters as follows: d = 2, 
r = 2, j s f = 0.1, E F = 2.45, 7sp G {0.05,0.1,0.2,0.25}. 
The figure shows a gate-induced spin current modula- 
tion whose voltage period, independent from 7 sp / s /, is 
given by r„ = ir/d. The same period is found in Fig llll 




e v a 



FIG. 10: (Color online) Spin currents J s vs eV g in gate config- 
uration 2 computed by fixing the model parameters as follows: 
d — 2, F = 2, 7 S / = 0.1, Ef = 2.45. The parameter 7„ p varies 
as 7s P G {0.05,0.1,0.2,0.25}. 
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where we analyze the spin current J s as a function of eV g 
fixing the model parameters as follows: d = 2, T = 2, 
7 sp = 0.1, E F = 2.45, 7s/ € {0.05,0.1,0.2,0.25}. In 
this case a peculiar form of spin-switch occurs. Indeed, 
at increasing values of jgf the curves show a change of 
sign of the spin current when the gate is tuned close to 
odd integer values. However, the possibility to obtain a 




-6 -5 -4 -3 -2 
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FIG. 11: (Color online) Spin currents J„ vs eV g in gate config- 
uration 2 computed by fixing the model parameters as follows: 
d = 2, F = 2, j S p = 0.1, Ef — 2.45. The parameter y s f varies 
as j sf e {0.05,0.1,0.2,0.25}. 

full sign reversal is determined by the constriction prop- 
erties (i.e. 7s/). To fully characterize the spin-switching 
phenomenon, in Fig |12l we report the spin current as a 
function of the constriction length d by fixing the model 
parameters as follows: eV g = —2.4, T = 2, ^ sp = 0.1, 
E F = 2.45, j sf G {0.05,0.1,0.2,0.25}. The figure shows 
a complicated beating-like pattern whose periods are de- 
termined by the ones in Eq. (|2"0")) . The interesting aspect 
in FigfT^] is that, in the parameters range investigated, 
there exist special values of d where the corresponding 
value of J s is independent on 7s j . This property implies 
that devices having a value of d close to these special 
points do not manifest the spin switching phenomenon. 

V. CONCLUSIONS 

We investigated the electrical switching of charge and 
spin transport in a topological insulator nanoconstriction 
in a four terminal device by means of a scattering field 
theory. The spin and charge switch of the edge chan- 
nels is caused by the coupling between edges states which 
overlap in the constriction and by the tunneling effects 
at the contacts of the quantum spin Hall bar and there- 
fore can be manipulated by tuning external applied gate 
voltages and by geometrical etching. We showed that 
the switching mechanism can be conveniently studied by 
electron interferometry involving the measurements of 
charge and spin currents in specific bias and gate con- 
figurations (the switching probability is analyzed in the 
Appendix [Bj . The device can operate in two gate config- 




FIG. 12: (Color online) Spin currents J s vs d in gate configu- 
ration 2 computed by fixing the model parameters as follows: 
eV g = -2.4, r = 2, 7s P = 0.1, E F = 2.45. The parameter 7s/ 
varies as j sf € {0.05,0.1,0.2,0.25}. 



urations defined as: (i) Vt = H = V g ; (ii) Vt = — V& = V g , 
with side gates Vt and V&. Furthermore, operating with 
the external dc bias Vi, the device can work in the fol- 
lowing configurations: (i) Charge-bias (CBC) defined as 
Vi = V 2 = V, V 3 = V 4 = 0; (ii) Spin-bias (SBC) defined 
as -Vi = V 2 = V,V 3 = V 4 = 0. The CBC and the SBC 
produce pure charge and pure spin current, respectively. 
Concerning the gate configuration (i), we showed that 
in CBC the Fabry-Perot interference maxima of J c vs 
eV g depend on the scattering phase S , while the quasi- 
periodic oscillating behavior is controlled by the coupling 
energy T. The above behavior is peculiar for a gapped 
spectrum of the quasi-particles, while for large enough 
transverse dimension W the maxima position does not 
depend on <p s , as in the case of massless Dirac parti- 
cles. A similar behavior is also observed in the J c vs d 
curves, where an oscillating dependence on the constric- 
tion length is present. The effect of cj) s on the transport 
properties is also evident in the SBC configuration. In 
particular, the J s vs eV g curves present an oscillating pat- 
tern reminiscent of the quantum interference, while the 
behavior of J s vs d curves is qualitatively similar to the 
one observed for the charge current in CBC. However, the 
spin current is strongly dependent on the spin-flipping 
tunneling value at the etching points x\^- In particular, 
the increasing of "f s f above a certain threshold value in- 
duces a change of sign of J s , the latter being important 
to characterize the interface properties of the device and 
for spintronics purposes. 

Concerning the gate configuration (ii), in the CBC we 
observe that the J c vs eV g curves do not depend on 
eV g . This behavior can be understood observing that the 
Fabry-Perot phase (f>FP responsible for the gate modula- 
tion of J c is zero. On the other hand, the J c vs d curves 
present a behavior similar to the one observed in the gate 
configuration (i). 

In SBC, the J s vs eV g curves present an oscillating behav- 
ior induced by an Aharonov-Bohm-like phase 4>ab 0, 
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which is strongly affected by the tunneling amplitudes 
7 sp and j s f. In particular for high values of j s f, a full 
electrical switching of the spin current can be obtained 
by tuning eV g in the interval [—3.9,-3.2] (see FigfTTj). 
Indeed, using the external gates, it is possible to tune 
the spin current from J s = —2 (eV g = —3.9) up to 
J s = 0.8 (eV g — —3.2), being this effect relevant for spin- 
tronics. Moreover, the J s vs d curves present a quantum 
beating-like structure caused by the presence of two (non- 
commensurate) periods of oscillation, namely t^ 2 ^ . Fi- 
nally, there exist special values of d for which the spin 
current is independent from "f s f (see Figfl2|). 
In conclusion, we have shown that the electrical switch- 
ing behavior offers an efficient and robust mechanism to 
manipulate topological edge states transport exploiting 
a four terminal set-up relevant for spintronics applica- 
tions. Furthermore, going beyond the technological im- 
plications, our analysis can be used to unveil the effects 
of a non-vanishing effective mass on the transport prop- 
erties of confined Dirac Fermions. 



Appendix A: Equation of motion and model 
assumptions 

Within the constriction [x\ < x < x 2 ), the equation 
of motion of ^(x) = {ipnf, 4>ri, i>LiY is governed by 
the Hamiltonian H = Hq + H c where for generality, we 
assume that apart from the term © a term with spin- 
flipping coupling of the form (j4|) with constant amplitude 
Tf is also present: 



ihd t ^(x) 



v + Tf r 

Tf v + r 

r V- -Tf 

r -r f V- 



(Al) 



with T>± — ^fihvFdx- Considering constrictions whose 
transverse dimension W allows edge states coupling, the 
functional form of the coupling 16 is T ss To exp(— XW) 
(see Fig fTS")) . while the flipping term Tf can be neglected 
for not too tight and not too long constrictions. As also 
evident from Ref.[13j, the spin precession induced by a 
very weak value of the spin-orbit interaction assumes a 
relevant role in the transport only for very long constric- 
tions characterized by L ps 900— 1900 nm which are quite 
long to maintain coherent transport. This comes from 
the fact that a slow spin precession needs a very long 
time (and thus a long distance) to completely change the 
spin polarization of a particle traveling along the sys- 
tem. In particular, a spin precession described by an 
angular frequency u>so requires a time r = it/luso to 
completely flip the electron spin, while the dwell time 
T~dw required to cover the constriction length L is given 
by L/vp. The above condition implies that the shortest 
constriction length to observe a complete spin flipping 
fulfills the relation L = ttvf/ujso- As a consequence, 
weak values of the spin orbit coupling (and thus small 
values of ujso) require longer constrictions. Motivated 
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FIG. 13: (Color online) Exponential dependence of F on W 
as deduced from Ref.[Tjj], the interpolation formula being V ~ 
F exp(-AW), with T = 7.75 meV and A = 0.018 nm" 1 . 



by these arguments, in this work we assume L < 1 /im 
and set Tf = along the constriction, which is equiva- 
lent to neglect the spin-orbit coupling, while we consider 
spin- flipping tunneling j s f at the etching points x = Xi/ 2 - 



Appendix B: Transition probabilities 

In this Appendix we show how topological edge states 
can be selectively switched in the nanoconstriction. In 
FigHH are shown the transition probabilities for an in- 
coming spin-polarized state at the upper edge (terminal 
2) to be reflected back to the lower edge (terminal 1), 
transmitted through the constriction in the same spin 
and edge state (terminal 3) or transmitted by swapping 
the edge and simultaneously flipping the spin (termi- 
nal 4). By using the parameters eV g — —2.45, T = 2, 
7s P = 0.1 (as given in FigjS]), and fixing the electrochem- 
ical potential at the bottom of the conduction band, one 
observes that transmission from a state up to a state 
down in terminal 4 is activated for longer nanoconstric- 
tion channels and higher values of the local spin flipping 
tunneling at the point contacts (see right down panel), 
while for lower values of j s f transmission along the same 
edge and spin state is favored for longer channels (see 
panel right up). For intermediate values of 7 S / both re- 
flection to the lower edge in the same spin state and trans- 
mission along the channel in terminal 3 and 4 takes place, 
the latter having almost equal probabilities (see middle 
panels). Since the spin-flipping tunneling is controlled 
by the local modification of the spin-orbit interaction by 
geometrical etching, our device can be used as spin tran- 
sistor with high fidelity. Let us note that differently from 
Ref.[l3|], here spin precession along the channel does not 
take place and is not relevant for this setup. 
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